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ABSTRACT 

We investigate the role of massive perturbers, such as stellar clusters or giant 
molecular clouds, in supplying low-angular momentum stars that pass very close to 
the central massive black hole (MBH) or fall into it. We show that massive perturbers 
can play an important role in supplying both binaries and single stars to the vicinity of 
the MBH. We discuss possible implications for the ejection of high velocity stars; for 
the capture of stars on tight orbits around the MBH; for the emission of gravitational 
waves from low-eccentricity inspiraling stars; and for the origin of the young main 
sequence B stars observed very near the Galactic MBH. Massive perturbers may also 
enhance the the growth rate of MBHs, and may accelerate the dynamical orbital decay 
of coalescing binary MBHs. 



1. Introduction 

There is compelling evidence that massive black holes (MBHs) lie in the centers of all galax- 
ies (e.g. (i 10;)), including the center of our Galaxy 

am. 

The MBH affects the dynamics and 
evolution of the galaxy's center as a whole and it also strongly affects individual stars or binaries 
that approach it. Such close encounters have been the focus of many studies and include a variety 
of processes such as destruction of stars by the MBH; capture and gradual inspiral of stars into 
the MBH, accompanied by the emission of gravitational waves (GWs); or dynamical exchange 
interactions in which incoming stars or binaries energetically eject a star tightly bound to the MBH 
and are captured in its place very near the MBH <£xl) - 

The interest in such processes is driven by their possible implications for the growth of MBHs, 
for the orbital decay of a MBH binary, for the detection of MBHs, for GW astronomy, as well as 
by observations of unusual stellar phenomena in our Galaxy, e.g. the puzzling young population 
of B-star very near the Galactic MBH , or the hyper- velocity B stars at the edge of the Galaxy (e.g. 
fiEE)), possibly ejected by 3-body interactions of a binaries with the MBH. 
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Here we focus on close encounters with the MBH whose ultimate outcome ("event") is the 
elimination of the incoming object from the system, whether on the short infall (dynamical) time, if 
the event is prompt (e.g. tidal disruption or 3-body exchange between a binary and the MBH), or on 
the longer inspiral time, if the event progresses via orbital decay (e.g. through GW emission). Such 
processes are effective only when the incoming object follows an almost zero angular momentum 
("loss-cone") orbit with periapse closer to the MBH than some small distance q. To reach the MBH, 
or to decay to a short period orbit, both the infall and inspiral times must be much shorter than the 
system's relaxation time t r . The fraction of stars initially on loss-cone orbits is very small and 
they are rapidly eliminated. Subsequently, the close encounter event rate is set by the dynamical 
processes that refill the loss-cone. 

The loss-cone formalism used for estimating the event rate (e.g.([z|)) usually assumes that the 
system is isolated and that the refilling process is 2-body relaxation. This typically leads to a low 
event rate, set by the long 2-body relaxation time. 

Two-body relaxation, which is inherent to stellar systems, ensures a minimal loss-cone re- 
filling rate. Other, more efficient but less general refilling mechanisms were also studied with the 
aim of explaining various open questions, or in the hope that they may lead to significantly higher 



event rates for close encounter processes (25). However, most of these mechanisms require special 
circumstances to work, or are short-lived. 

Here we explore another possibility, which is more likely to apply generally: accelerated re- 
laxation and enhanced rates of close encounters driven by massive perturbers (MPs). Efficient 



relaxation by MPs were first suggested in this context by Zhao, Haehnelt & Rees (1391) as a mech 



anism for establishing the M./a relation d 1 Oh by fast accretion of stars and dark matter. Zhao 
et al. also noted the possibility of increased tidal disruption flares and accelerated MBH binary 
coalescence due to MPs. In this study we investigate in detail the dynamical implications of relax- 
ation by MPs. We evaluate its effects on the different modes of close interactions with the MBH, 
in particular 3-body exchanges, which were not considered by Zhao et al. and apply our results 
to the Galactic Center (GC), where observations indicate that dynamical relaxation is very likely 
dominated by MPs. 



2. Loss-cone refilling 

In addition to stars, galaxies contain dense massive structures such as molecular clouds, open 
clusters and globular clusters with masses up to 1O 4 -1O 7 M . Such structures can perturb stellar 
orbits around the MBH much faster than 2-body stellar relaxation (hereafter "stellar relaxation"), 
provided they are numerous enough. The minimal impact parameter still consistent with a small 
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angle deflection in the MP-star scattering is 6 min = GM p /v 2 (the capture radius), where v is of the 
order of the local velocity dispersion a. The total rate of scattering stars into the loss cone, T, is 
obtained by integrating dT/dM p db over all MP masses and over all impact parameters between 
r. The relaxation rate due to all MPs at all impact parameters is then 
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where log A = log(6 max /6 m in) is the Coulomb logarithm (here the dependence of log A and v 
on M p is assumed to be negligible). This formulation of the relaxation time is equivalent to its 
conventional definition (Q) as the time for a change of order unity in v 2 by diffusion in phase space 
due to scattering, t r ~ v 2 /D(v 2 ), where D(v 2 ) is the diffusion coefficient. If the stars and MPs 
have distinct mass scales with typical number densities A 7 * and N p and rms masses {M 2 } 1/2 and 
(M 2 ) 1/2 ((M 2 ) = J M 2 (dN/dM)dM/N), then MPs dominate if the ratio of the 2nd moments of 
the MP and star mass distributions, fi 2 = N p (M 2 ) / A* (M 2 ), satisfies /i2^> 1. 

The central ~ 100 pc of the GC contain 10 s — 10 9 solar masses in stars, and about 10 6 — 10 8 
solar masses in MPs such as open clusters and GMCs of masses io 3 - io 7 M H;0;HH;li). 
An order of magnitude estimate indicates that MPs in the GC can reduce the relaxation time by 
several orders of magnitude, 
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This estimate is borne by more detailed calculations (Fig. [Hand tabled)), using the formal defi- 
nition t r = v 2 /D(v 2 ) with M±p± -> / (dN p /dM p )M?dM p (Q). In our calculations we follow the 
Fokker-Planck approach to the loss-cone problem (17|), where we recalculate the diffusion coeffi- 
cients obtained from stellar two body relaxation (Q) by taking into account the contribution from 
the much more massive MPs (for details see d^)). 

Although some of the assumptions concerning the loss cone formalism are not necessarily 
valid in the case of MPs, the loss-cone formalism can be generalized to deal with MPs in an ap- 
proximate manner with only few modifications (|32|) . Thus to a good approximation our calculation 
follows the usual loss cone treatment (e.g. (|24l) ) where the mass of the stars is replaced by that of 
the MPs, and the integration over the energies of deflected stars (mapped to their distance from the 
MBH) is done only for regions where MPs exist. 



3. Massive perturbers in the Galactic Center 



MPs can dominate relaxation only when they are massive enough to compensate for their 
small space densities. Here we consider only MPs with masses M p > 10 3 M Q . Such MPs could 
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Fig. 1. — Relaxation time as function of distance 
from the MBH, for stars (solid line) and for each of 
the 2 MP models separately, as listed in table Q] mas- 
sive GMCs (dashed-dotted line) and intermediate GMCs 
(dashed line). The discontinuities are artifacts of the as- 
sumed sharp spatial cutoffs on the MP distributions. 2- 
body stellar processes dominate close to the MBH, where 
no MPs are observed to exist. However, at larger dis- 
tances massive clumps (at 1.5 < r < 5 pes) and GMCs 
(at 5 < r < 100) are much more important. 



be open or globular stellar clusters or molecular clouds of different masses, in particular giant 
molecular clouds (GMCs). Observations of the Galaxy reveal enough MPs to dominate relaxation 
in the central 100 pc. We adopt here a conservative approach, and include in our modeling only 
those MPs that are directly observed in the Galaxy. Based on the observations of such MPs we 
devise several possible models for the MPs in the GC, detailed in tabled! As these observations (0; 

show that the MPs population is dominated by the GMCs, we consider only GMCs 



and gaseous clumps in our models. The observed MP species vary in their spatial distributions 
and mass functions, which are not smooth or regular. For our numeric calculations, we construct 
several simplified MP models (table [B that are broadly based on the observed properties of the 
MPs, and assume that the spatial distribution of MPs follows that of the stars (apart of the inner 
radius cutoff where no MPs are observed). 

The three MP models in table Q] Stars, GMC1 and GMC2, represent respectively the case of 
relaxation by stars only, by heavy GMCs and light GMCs. 



Table 1: Massive perturber models 



Model 


r (pc) a 


N p 


M p (M ) 


p 


R p (pc) 


H2 b 


GMC1 


5-100 


100 


10 4 -10 7 


1.6 


5 


3xl0 5 




1.5-5 


30 


10 3 -10 5 


0.9 




6000 


GMC2 


5-100 


100 


10 3 -10 6 


1.6 


5 


3xl0 3 




1.5-5 


30 


10 2 -10 4 


0.9 




60 


Stars 


5-100 


2xl0 8 


1 




~0 


1 


Stars 


1.5-5 


6xl0 6 


1 




~0 


1 


a N p (r)ocr 2 assumed. 

b b /i 2 = N p {Ml) /N* (M?) , where <M 2 ) =/ M 2 (dN/dM)dM/N. 
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4. Massive perturber-driven interactions with a MBH 

The maximal differential loss-cone refilling rate, which is also the close encounters event rate, 
dT/dE, is reached when relaxation is efficient enough to completely refill the loss cone during one 
orbit. Further decrease in the relaxation time does not affect the event rate at that energy. MPs can 
therefore increase the differential event rate over that predicted by stellar relaxation, only at high 
enough energies, E > E c (equivalently, small enough typical radii, r < r c ), the critical energy 
(radius), separating the full and empty loss cone regimes, where slow stellar relaxation fails to 
refill the empty loss-cone. The extent of the empty loss-cone region increases with the maximal 
periapse q, which in turn depends on the close encounter process of interest. For example, the 
tidal disruption of an object of mass M and size R occurs when q < r t , the tidal disruption radius, 

1 /3 

rt — R (M./M) ' . This approximate disruption criterion applies both for single stars (M = M*, 
R = R+) and for binaries, where M is the combined mass of the binary components and R is the 
binary's semi-major axis, a. Stellar radii are usually much smaller than typical binary separations, 
but stellar masses are only ~ 2 times smaller than binary masses. Binaries are therefore disrupted 
on larger scales than single stars. In the GC this translates to an empty (stellar relaxation) loss-cone 
region extending out to r s c ~ 3 pc for single stars and out to r h c > 100 pc for binaries. In the GC 
r MP ^ r c ^ r c (where tmp is the smallest distance from the MBH where MPs are observed) , and 
so MPs are expected to increase the binary disruption rate by orders of magnitude, but increase the 
single star disruption rate only by a small factor. Since for most Galactic MP types r M p > r h (the 
radius of influence of the MBH), the disruption rate is dominated by stars near r M p. For example, 
when the loss -cone is empty, ~ 50% of the total rate is due to MPs at r < 2rMp; when the loss -cone 



is full, ~ 75% of the total rate is due to MPs at r < 2r M p (see details in (|32l) ) 



4.1. Interactions with single stars 

Clusters, GMCs and gas clumps in the GC are abundant only beyond the central r MP ~ 1.5 pc, 
whereas the empty loss-cone regime for tidal disruption of single stars extends only out to r s c ~ 3 



pc. For inspiral processes such as GW emission, r c is ~ 100 times smaller still (|19l) . The effect of 
such MPs on close encounter events involving single stars is thus suppressed (weaker tidal effects 
by MPs at r > r s c are not considered here). This is contrary to the suggestion of Zhao et al. (|39l) . 
who assumed that the effect of MPs fully extends to the empty loss-cone regime. We find that the 
enhancement of MPs over stellar relaxation to the single stars disruption rate is small, less than a 
factor of 3, and is due to stars scattered by gas clumps in the small empty-loss cone region between 
tmp ~ 1.5 pc and r s c ~ 3pc. A possible exception to this conclusion is the possible existence of 
IMBHs population, not modeled here. 



-6- 



4.2. Interactions with stellar binaries 

The empty loss cone region for binary-MBH interactions extends out to > 100 pc because of 
their large tidal radius. On these large scales MPs are abundant enough to dominate the relaxation 



processes. Here we focus on 3-body exchange interactions (1171:1381). which lead to the disruption of 
the binary, the energetic ejection of one star, and the capture of the other star on a close orbit around 
the MBH. The event rate is highly dependent on the unknown binary fraction in these regions. 

The binary fraction and typical binary semi-major axis depend on the binary mass, and on 
the rate at which binaries evaporate by encounters with other stars. This depends in turn on the 
stellar densities and velocities, and therefore on the distance from the MBH. We take these factors 
into account and estimate in detail the 3-body exchange rate for MP-driven relaxation. The rate 
is proportional to the binary fraction in the population, which is the product of the poorly-known 
binary IMF in the GC and the survival probability against binary evaporation. 

The capture probability and the semi-major axis distribution of captured stars were estimated 



by simulations (|18l : I38|) . Numeric experiments indicate that between 0.5-1.0 of the binaries that 
approach the MBH within the tidal radius r t (a) are disrupted. Here we adopt a disruption efficiency 
of 0.75. The harmonic mean semi-major axis for 3-body exchanges with equal mass binaries was 



found to be (|18|) 



, v ( M. \ 2/3 / M. \ 1/3 

where a is the semi-major axis of the infalling binary and ai that of the captured star (the MBH-star 
"binary"). Most values of a\ fall within a factor 2 of the mean. This relation maps the semi-major 
axis distribution of the infalling binaries to that of the captured stars: the harder the binaries, the 
more tightly bound the captured stars. The periapse of the captured star is at r t , and therefore its 



eccentricity is very high (|18c |27|) . e = 1 — r t /ai — 1 — 1.8(M bin /M.) 1 / 3 > 0.98 for values typical 
of the GC. 

We now consider the implications of 3-body exchange interactions of the MBH with old 
(t* >tff) binaries and massive young (£* < 5 x 10 7 yr) binaries. 

The properties of binaries in the inner GC are at present poorly determined. We use the 
period distribution of Solar neighborhood binaries for old low mass binaries dsl) and young massive 



binaries (|21|) . The total binary fraction of these binaries is estimated at / bin ~ 0.3 for low mass 
(|22l) binaries, and f hin ~ 0.75 for massive binaries (|2ll) . Adopting these values for the GC, the 
total binary disruption rate by the MBH can then be calculated by integrating diVbin/da over the 
binary a distribution and over the power-law stellar density distribution of the GC from the minimal 
radius where such binaries exist (for young binaries we assume an inner cutoff at 1 .5 pes, where 



such young stars are not observed (|3ll) ). up to 100 pc (|13l) . Table © lists the capture rates for the 
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different perturber models, assuming a typical old equal-mass binary of M^ m = 2 M Q , or young 
equal-mass binary of M bin = 15 M . 

The old, low-mass binary disruption rate we derive for stellar relaxation alone is ~ 5 x 



10 yr , ~ 5 times lower, but still in broad agreement with the result of Yu & Tremaine (38). 
Their rate is somewhat higher because they assumed a constant binary fraction and a constant 
semi-major axis for all binaries, even close to the MBH, where these assumptions no longer hold. 

MPs increase the binary disruption and high- velocity star ejection rates by factors of ~ 10 1 " 3 
and effectively accelerate stellar migration to the center. Thus we expect an increase in the number 
of stars captured close to the MBH and consequently a higher event rate of single star processes 
such as tidal disruption, tidal heating and GW emission from compact objects, in particular from 



compact objects on zero-eccentricity orbits (I27|) . MPs may be implicated in the puzzling presence 
of a cluster of main sequence B-stars (4 < M* < 15 M ) in the inner ~ l" (~ 0.04 pc) of the 
GC. This so-called "S-cluster" is spatially, kinematically and spectroscopically distinct from the 
young, more massive stars observed farther out, on the ~ 0.05-0.5 pc scale, which are thought to 



have formed from the gravitational fragmentation of one or two gas disks (|31l) . There is however 
still no satisfactory explanation for the existence of the seemingly normal, young massive main 
sequence stars of the S-cluster, so close to a MBH (see review of proposed models by Alexander 



1); also a recent model by Levin (|23l) . and in these proceedings). 



Here we revisit an idea proposed by Gould & Quillen (15), that the S-stars were captured 
near the MBH by 3-body exchange interactions with infalling massive binaries. Originally, this 
exchange scenario lacked a plausible source for the massive binaries. We suggest that MP-driven 
3-body exchanges can serve as a source for binaries, as they increase the rate of infall of young 
field binaries to the MBH. Such young field binaries should exist in the GC, taking into account 
the ongoing star formation in the central ~ 100 pc, which implies the presence of a large reservoir 
of massive stars there. Such stars are indeed observed in the central few x 10 pc both in dense 



Table 2: Total binary disruption rate and number of captured young stars 



Model 


Disruption rate (yr x ) 


Young Stars" 


Young Stars 




r<0.04pc r<0.4pc 


r <0.04pc 


0.04 < r<0.4pc 


GMC1 


1 x 10" 4 2.8 x 10" 4 


33.5 


3.8 


GMC2 


2.1 x 10~ 5 6.4 x 10~ 5 


5 


0.24 


Stars 


3.4 x 10~ 7 5.3 x 10~ 7 


0.15 


0.003 


Observed 


? ? 


10 - 35 b 


? 



"Main sequence B stars with lifespan t < 5 x 10 7 yr . 

b ~ 10 stars with derived a< 0.04 pc. >30 stars are observed in the area. 
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clusters and in the field |29|) . It is plausible that a high fraction of them are in binaries. 



We assume star formation at a constant rate for 10 Gyr with a Miller-Scalo IMF (1261) . and use 



a stellar population synthesis code (1351) with the Geneva stellar evolution tracks (1341) to estimate 
that the present day number fraction of stars in the S-star mass range is 3.5 x 10~ 4 (and less than 
0.01 of that for M* > 15 M stars). Note that if star formation in the GC is biased toward massive 



stars (1361) . this estimate should be revised upward. We adopt the observed Solar neighborhood 



distribution of the semi-major axis of massive binaries, which shows that massive binaries are thus 



typically harder than low-mass binaries (2jJ), and will be tidally disrupted closer to the MBH and 
leave a more tightly bound captured star. 

We represent the massive binaries by one with equal mass stars in the mid-range of the S-stars 
masses, with M hin = 15 M Q and i*(7.5 M ) ~ 5 x 10 7 yr, and integrate over the stellar distribution 
and the binary a distribution as before, to obtain the rate of binary disruptions, T, the mean number 
of captured massive stars in steady state, TV* = Tt+, and their semi-major axis distribution (Eq. 
[3]). Table © compares the number of captured young stars in steady state, for the different MP 
models, on the r < 0.04 pc scale (the S-cluster) and 0.04 < r < 0.4 pc scale (the stellar rings) with 
current observations (^liD). 

The number of captured massive stars falls rapidly beyond 0.04 pc (table [2]) where the S- 
cluster is observed because wide massive binaries are rare. This capture model thus provides a 
natural explanation for the central concentration of the S-cluster (Fig [2]). The absence of more 
massive stars in the S-cluster (M* > 15 M , spectral type O V) is a statistical reflection of their 
much smaller fraction in the binary population. Figure © and table © compare the cumulative 
semi-major axis distribution of captured B-stars, as predicted by the different MP models, with the 



total number of young stars observed in the inner 0.04 pc (~ 35 stars (|9|; 114J; 1311) . Of these, only 
~ 10 have full orbital solutions (in particular a and e) at present. The numbers predicted by the 
MP models are consistent with the observations, unlike the stellar relaxation model that falls short 
by two orders of magnitude. 

The binary capture model predicts that captured stars have very high initial eccentricities. 
Most of the solved S-star orbits do have e > 0.9, but a couple have e ~ 0.3-0.4 (@). Normal stellar 
relaxation is too slow to explain the decrease in the eccentricity of these stars over their relatively 



short lifetimes. However, the much faster process of resonant relaxation (1331) may be efficient 



enough to randomize the eccentricity of a fraction of the stars, and could thus possibly explain the 



much larger observed spread in eccentricities (|20h . 



The companions of these captured stars are ejected from the GC at high velocities 

(0). Con- 
sequently, our model predicts the number of young hyper velocity stars ejected from the GC (such 
as observed (0;0;S)) to be comparable to that of observed B-stars in the S-cluster, in agreement 
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with current observation based estimations 
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Fig. 2. — Cumulative number of young B-stars 
in the GC as predicted by the MP models and 
by stellar two-body relaxation (listed in table 
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served. 



5. Summary 



We presented here the results of a study of the effect of massive perturbers near the MBH in 
the GC. We have shown that current observations of MPs such as GMCs and clusters indicate that 
they dominate relaxation process in the GC, where they exist, and they are much more important 
than stellar two-body relaxation processes by stars which are usually considered. We have used 
the loss cone formalism in order to analyze the importance of this mechanism in generating low 
angular momentum stars and binaries to interact with the MBH. We have have computed the rates 
of these interactions, and showed that they can be highly important for the process of binaries 
disruption by the MBH and its consequent outcomes. The origin of the young massive B-stars at 
the central arcsecond of the GC can be explained by the capture of stars from binaries disrupted by 
the MBH after being scattered by MPs. Our calculations also show that some of the companions of 
these captured stars could be ejected at high velocities, thus explaining the observations of young 
massive high velocity stars observed recently. In addition MPs increase the number of compact 
stars captured close to the MBH, and thus increase the emission of zero-eccentricity GWs. We 
also suggest that such MPs may also help solve the last parsec problem of coalescing BMBHs. 
Although we focused on the GC of the milky-way, our results can be easily extended to MPs near 
other MBHs. 
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